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^ , Distances play important roles in cosmological observations, especially in gravitational 



lens systems, but there is a problem in determining distances because they are defined 



' in terms of light propagation, which is influenced gravitationally by the inhomogeneities 

in the universe. In this paper we first give the basic optical relations and the definitions 
of different distances in inhomogeneous universes. Next we show how the observational 
relations depend quantitatively on the distances. Finally, we give results for the frequency 
distribution of different distances and the shear effect on distances obtained using various 
methods of numerical simulation. 



■ SI. Introduction 

I In optical relations among observed quantities, distances such as the luminos- 

0^ . ity distance and the angular diameter distances play an important role. They are 

clearly defined in the homogeneous Priedmann-Lemaitre-Robertson- Walker model 
I (Weinberg, EI Schneider et aLW) owing to the simple nature of light propagation in 

this case. In inhomogeneous universes, however, their behavior is complicated, due 
to gravitational lens effect which implies that light rays are deflected gravitationally 
^1 by an inhomogeneous matter distribution. On the other hand, we also use distances 

■ to interpret the structure of gravitationally lensed systems. 
To correctly treat distances in inhomogeneous universes, it is necessary first to 

^ have a reasonable formulation for the dynamics describing local matter motion and 

■ optics and clarify the validity condition of the formulation. A set of fluid dynamical 
equations and the Poisson equation in the cosmological Newtonian approximation 
was introduced and discussed by Nariaiu^ and IrvineD* under the conditions 

\^\ < 1, {v/cf < 1, L/Lh < 1, (M) 

where v, L and Lh are the Newtonian gravitational potential, matter velocity, 
the characteristic size of inhomogeneities and the horizon size ~ ct, respectively, and 
the spacetime is expressed as 

ds^ = -(1 + 2^)c^dt^ + (1 - 2^)a^{t)[dx^ + (j^{x)dn^], (1-2) 

where a{t) is the scale-factor, a{x) = sinx,XiSinhx for the background curvature 



*•* E-mail address: tomita@yukawa.kyoto-u.ac.jp 
**-* E-mail address: asada@phys.hirosaki-u.ac.jp 
**•) E-mail address: hamana@astr.tohoku.ac.jp 



156 



K. Tomita, H. Asada and T. Hamana 



k = 1,0,-1, respectively, and dO^ = dO^ + sit? O^p^ . The above fluid dynamical 
equations can describe the nonlinear local motion, while the gravitational field is 
linear with respect to ^. The extension of the above cosmological Newtonian treat- 
ment to a post-Newtonian treatment was performed by Futamase,^ TomitaS and 
Shibata and Asada. Futamase showed that the condition 

e^/K<^l (e^ ~ ^ and K ~ L/Li^) (1-3) 

is necessary for the higher-order expansion to be possible and formulated the spa- 
tial averaging and the back-reaction to the background. Moreover, Futamase and 
Sasaki eP investigated the validity of light propagation in the cosmological Newtonian 
iterative approximation and discussed the distance problem. 

In an empty region as the limiting inhomogeneous case, distances exhibit be- 
havior very different from those in homogeneous models (the Friedmann distances). 
In the special case without tidal shear from surrounding regions, the so-called Dyer- 

toeder angular diameter distance was derived by Zerdovich,l2P Dashevskii and Slysh 
\ and Dyer and Roeder.0) In a non-empty region with a constant matter den- 
sity pm but no tidal shear, we have the generalized Dyer-Roeder distance with the 
dumpiness (or smoothing) parameter a, which is defined as 

Pml Py = a = const (1'4) 

for the Friedmann density pp. The observational results derived from the optical 
relations depend on whether we use the Friedmann distances or the Dyer-Roeder 
distance. Quantitative estimates for these differenceand the effect of the cosmolog- 
ical constant have been studied by Fukugita et al.ll^) and Asada. On the other 
hand, it is important to determine what distances are most applicable and what 
value of the above parameter a is best, in realistic inhomogeneous models. Kasai 
et al.© and Watanabe and Tomita© numerically calculated the frequency distri- 
bution for generalized distances in simple models in which particles are distributed 
randomly. Recently, Tomita derived this distribution in more realistic inhomoge- 
neous models generated using the A^-body simulation with the CDM spectrum. The 
general result is that the average value of a is nearly 1 and its dispersion decreases 
with the increase of the redshift z, though it is ~ 1 for z = 0.5. 

Another interesting topic is that involving the role of the shear term and the 
Ricci and Weyl focusing terms (in the optical scalar equation 0^) in the 
distances. To this time, the shear effect has been discussed by Weinberg, 
abe et al.,c3) Watanabe and Sasaki, EiP and Nakamura,E3) and its focusing effect 
has recently been studied numerically by Hamana using a Monte Carlo simulation, 
taking into account small-scale inhomogeneities. 

In this review paper, basic optical relations and the definition of distances are 
first given in §2, the lensing relations in the generalized Dyer-Roeder distance are 
derived in §3 (by Asada), the statistical behavior of distances analyzed in numerical 
simulation is described in §4 (by Tomita), and the shear and focusing effects are 
discussed in §5 (by Hamana). 



iehavior of 
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Distances in Inhomogeneous Cosmological Models 



157 



§2. Optical scalars and the definition of distances 

2.1. Geometry of ray bundles 

Rays are expressed as = x^{y^,v), where v is an affine parameter. For each 
ray, the have constant values (C*). The wave vector k^^ = dx^ jdv satisfies the 
nuh condition and null-geodesic equation 

}eko, = 0, k'',/3k'^ = 0. (2-1) 

For two rays with = C* and = C* + 6C^, the connection vector is 

fe^ = {dx''/dy')6C\ (2-2) 

If we define a dot differentiation by {m")' = m^.f^k^ for an arbitrary vector m^, we 
obtain from Eq. ( |2-2D 

= A;";^fe^, (2-3) 

and then (Jordan et al., ii) Sachs©) 

{kaSx^y = 0. (2-4) 

Now let us consider the situation in which on a screen (at a point Pq) an observer 
sees the shadow formed by a source object (at a point Pg). The connection vector 
6±x°' vertical to it" and k" at the point Pg is 

S^x'^ = h'^dx'^, (2-5) 

where 



k"kf3 k"u(3 + u^kp 

The hp satisfy the relation 



h'^h^ = h'^, hlu'^ = h'^k^ = and /i^ = 2. (2-7) 

If 6x" is the vector obtained by parallel-transporting d^x" from Pg along the ray, 
we have 

= d^x'^ + J {d^x'^ydv, (2-8) 
and the connection vector 5±x" vertical to u" and /c" at the point Pq is 

fo" = h'^Sx^ = (5^a;" + J h'^{6±x'^ydv. (2-9) 

The length 51 and the angle ai2 are expressed as follows using the connection vector 
in the plane vertical to /c" and the screen velocity : 

61 = {go^p~5x"~5x^f/\ cosai2 = [{5x'^)i{5x%]/[{5l)i{5l)2]. (2-10) 

Here u° is parallel-transformed along the ray from at Pg. 
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2.2. Optical scalar s 

Using Eqs. (|2-3| ) and (^) we obtain 

where 

The tensor A^p can be uniquely split as 



(2-11) 
(2-12) 



00 



a 



1 



2 

a/3 



(a/3) 



(2-13) 



where 0, a and w are optical scalars representing the expansion, shear and rotation, 
respectively, of ray bundles. In geometric optics which we assume in the following, 
the rotation vanishes, because is a gradient vector. By the transformation oiv, 
9, a and to transform, but 6dv, adv and Lodv are invariant (Jordan et al.,© SachsEl^). 
The evolution equations for 6 and a are 



dv ~ l'^-"" 

da 

dv 



-C - 2ea, 



(2-14) 



(2-15) 



where TZ = R^i^k^k'^, and C is expressed in terms of the Weyl tensor as C = 
Coii3'^\k'^k"'t^t^ , with a complex null vector satisfying t°'ta = k'^ta = 0,t'^ta = 1- 
As can be seen from Eqs. ( p-14D and ( 2-15 ) there two terms causing the focusing of 
ray bundles. One is the Ricci focusing term TZ, proportional to the matter density, 
and the other is the Weyl focusing term C, connected with the shear. 

2.3. Definition of distances 

The length of the shadow of the interval between two rays in the observer plane 
is given by 

{dlf = gaphx^hx''. (2-16) 

Then we have 

d{6l) 



51 



A^pe'^e'^dv = {9 + (T^pe'^e'^)dv, 



(2-17) 



where e° = 6±x'^/6l and ^a/je^e^ = 1. Next, let us consider the area of the cross- 
section of a ray bundle given by 



2tt 



{Slfda. 



(2-18) 
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If the deviation of the cross-section from a circle is small, we obtain 

d{5A) = {5Aydv = dv5A [^''^^da (2-19) 

Jo ol 

or 

^ = 2edv, (2-20) 

because Jf^^ Gape'^e^da vanishes. 

Here we define two kinds of angular diameter distances (the linear angular di- 
ameter distance Dia and the area angular diameter distance -DaA) proportional to 
51 and (5^4)^/^, respectively, as 



dv 

dDg^A 



+ a„;3e"e^) Aa, (2-21) 
0D,A, (2-22) 



dv 

where we consider ray bundles with 6 = oo at the observer point Pq- In the case 
of no shear, Dia and DaA are equal, but generally they are different. Their average 
values are equal if the term aafse'^e^ is cancelled out in the averaging process. In the 
situation that = oo at the source's point Pg, we obtain the luminosity distance Dl 
satisfying 

^=0D^. (2-23) 

dv 

The relation between D^a and I?l is proved to be 

Dl = {1 + zfD^A (2-24) 
by Etheringen, where z is the redshift given by the relation 1 + z = (u^fca) source 

/{u ^a) observer- 

The solutions of Eqs. ( ^-21 ) and (2-22) are generally complicated, but they can 
easily be obtained in the special case in which (1) the density is spatially constant, 
(2) there is no shear, and (3) the affine parameter in the Friedman background can 



be used. From Eqs. (|2T4| ) and (|2-22| ) we obtain 

'^-\nO., (2.2.) 

where Da = Aa = Ds^A, 7^ = SnGap/a^it) = 3(1 + z)^af2o, i7o is the total density 
parameter, and a is the smoothing (dumpiness) parameter. The affine parameter v 
is related to z as 

- = (1 + zf[{l + f2oz){l + zf - Xoz{2 + (2-26) 
dv 

where Aq is the normalized cosmological constant. The boundary condition for Da 
at epoch z = zi is 

„ dDA / , , , X c aiz-i ) 

Da = 0, — ^ = c{dt dz),=,, = -TT^^- 2-27 
az Hq ao 
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The solution of Eq. ( p-26| ) is called the Friedmann distance for a = 1, the Dyer- 
Roeder distance for a = 0, and the generalized Dyer-Roeder distance for arbitrary 
a. 

For the analyses of cosmological lens systems, we often use the lens equation 
;3 = 0_4G%^4 /rf'^'^(^')ln^^, (2-28) 



Dos 06 

where (3 and 6 are the angular position vectors (as seen by the observer) of the 
image and source, respectively, relative to the lens, 2J{9') is the surface mass density 
of the lens on the lens plane, 9c is an arbitrary constant angle, and Dol,Dos and 
Dls are the angular distances between the lens and the observer, the source and 
the observer, and the lens and the source, respectively. This equation has so far 
been derived only from intuitive geometrical considerations with use of the thin lens 
approximation, but it is not clear what distances should be used. In order tcLderive 
the lens equation from cosmological equations of light propagation, Sasaki E3) used 
the equation of geodesic deviation. It is obtained from Eq. ( |2-3D for the connection 
vector: 

(fe^)- = -R'^^^p^k'^dx'^. (2-29) 

He solved this equation in the case of an ideal light path, which is separated into 
three regions: the homogeneous, shearless region I (between the observer and the 
lens object), the region II including the lens object, and the homogeneous, shearless 
region III (between the lens object and the source). The light path in regions I and 
III is expressed by the generalized Dyer-Roeder distance, and the deflection of light 
rays in region III is determined using the thin lens approximation. The resulting lens 
equation reduces to the usual one with the generalized Dyer-Roeder distances. In the 
case that in the regions I and III there are inhomogeneous matter distributions and 
the shear effect is not negligible, however^the usual expression of the lens equation 
cannot be used, as was shown by Sasaki. ^•'^ 

In the multi-lens-plane method,!)^-!^ inhomogeneities as lens objects are 
assumed to be only in the lens planes, and hence the use of the lens equation in 
this method is consistent with the above assumption that the regions I and III are 
homogeneous and shearless. However the neglection of gravitational forces due to 
the difference of the projected matter distribution from the real distribution may be 
comparable with the neglection of weak forces from distant matter distribution. 



§3. Distances and lensing relations 

3.1. Observation and distances in gravitational lensing 

There are soine methods to determine cosmological parameters by using grav- 
itational lenses. El*'@^'©'©'il''i3)'@'li* Most of them concern the following three 
typical observational quantities: (1) the bending angle, (2) the lensing statistics and 
(3) the time delay. It is of great importance to clarify the determination of cosmolog- 
ical parameters through their observations in the realistic universe. In particular, it 
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has befiiL discuaseithBt Jmhm^ the universe may affect tlie cosmological 

testsB'i)'0)'0)©'E3)'Ei'i§M'B 

In this section, we use the so-called Dyer-Roedsr angular diameter distance in 
order to take account of the inhomogeneities.lill''E3)'Ei) We can consider this distance 
in two different cases, that of the so-called filled beam in the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe, and that of the so-called empty beam, when the 
right ray propagates through the empty region. For comparison with the filled beam, 
the empty beam has been freqiiently used and studied numerically in the literature 
(for instance, Fukugita et al.ta^'EJ)). However, it has not been clarified whether 
the observed quantities and/or the cosmological parameters for the arbitrary case 
of the dumpiness parameter are bounded between those for the filled beam and the 
empty beaiiL Moreover, numerical investigations have fixed redshifts of the lens and 
the source, though the effect of the dumpiness on the observable depends on 
the redshifts of the lens and the source. Therefore, it is important to clarify how 
the observation of gravitational lensing depends on all the parameters (the density 
parameter, cosmological constant, dumpiness parameter and redshifts of the lens and 
the source). For this reason, we derive the dependence on these parameters. 

3.2. Distance combinations in gravitational lenses 

(1) bending angle 

The lens equation is written asH^ 

(3 = e--^cx. (3-1) 

Here, /3 and are the angular position vectors of the source and image, respectively, 
and a is the vector representing the deflection angle. The effective bending angle 
{Di^Q,/ Dqi^)ol appears when we discuss the observations concemmg the angle such 
as the image separation and the location of the critical line. £3 'B' Hence the ratio 
Dls/Dqs plays an important role in the discussion of observations concerning the 
angle. It has been argued that, in calculating the bending angle, the density along 
the line of sight should be subtracted from the density of the lens object. l3* However, 
we assume that the density of the lens is much larger than that along the line of 
sight, so that the effect of the dumpiness on a can be ignored. Thus, we consider 
only the ratio -DlsZ-Cqs in the following. 

(2) lensing statistics 

The differential probability of lensing events is^'i^ 

dr = ani^dl, (3-2) 

where riL is the number density of the lens, dl is the physical length of the depth 
and a is the cross section, proportional to Dqi^Di^s/ Dqs- Since dl depends only on 
the cosmological parameters in the FLRW universe, we investigate the combination 
DolDi^s/ Dqq, in order to take account of the dumpiness of the matter. 

(3) time delay 
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The time delay between two images A and B is written as©'!) 

At^B = ' + '^''0^0^ f\e.( ^^±^ - c.{e)) , (3-3) 



where qa and qb are the bending angles at the images A and B, respectively. 

3.3. Monotonia properties 

It is assumed that the affine oarameter in the Dyer-Roeder distance is the same as 
that in the FLRW universe, 0^'© namely Eq. ( |2-26| ). Since a represents the strength 
of the Ricci focusing along the line of sight, Jihe DR angular diameter distance is a 
decreasing function of a for a fixed redshift,t3' that is to say, 

-DoL(ai) > -DoL(a2) for ai < 02- (3-4) 

(1) Dls/^os ^ 
It has been shown that the distance ratio -Dls/-^^os satisfies £3* 

^^("1) < ^^^(02) for ai < 02. (3-5) 
Dos Uos 

This is shown as follows. For fixed zs, Qq and Aq, the ratio i^Ls/^os can be 
considered as a function of zl, Xa{zi). We define Ya{zi) as -DslZ-DoSj where -DsL-is 
the Dyer-Roeder distance from the source to the lens. Owing to the reciprocity, EIP 
we obtain 

^a(^L) = \^X^{ZI^. (3-6) 

Since depends on zl only through DsL, it obeys the equation 

ra(2L) + |(l + ^L)'af2o^«(zL) =0, (3-7) 



where is an affine parameter at the lens. Let us define the Wronskian as 

WiY^^,Y„) = Yj^-Yj^. (3.8) 

Then, we obtain 

4-WiYo,.,Yao) <^ for ai<a2. (3-9) 
Since both Ya^ and Ya^ vanish at zl = zq,, we obtain 

w^(y,,(zs),y,,(zs)) = 0. (3-10) 



From Eqs. ( |3-9| ) and ( 3-10 ), we find 

W{Yo^,,Y^,)>Q, (3-11) 

where we used the fact that the affine parameter v defined by Eq. ( 2-26[ ) is an 
increasing function of z. Equation ( |3T1| ) can be rewritten as 



Distances in Inhomogeneous Cosmological Models 



163 



Since always becomes 1 + the observer, we find 

5^«2(^L = 0) 



From Eqs. ( |34^ ) and ( Fl3| ), we obtain 

^ > 1. (3-14) 

Thus, from Eq. (|3-6| ), Eq. (|3-5| ) is proved. 

From Eqs. ( |3-lD and ( |3-5D , we see the image separation as well as the effective 
bending angle increases with a. 
(2) DolDls/Dos 

Next let us prove that Dql^ls/^os increases monotonically with a. We fix i7, 
A, zl and zs- Then it is crucial to note that the distance from the lens to the soiirce 
can be expressed in terms of the distance function from the observer, D{z), asH 



Dls = Tr(l + ^i^)DolDos / (3-15) 
Ho Jdl D{zy 



dv 

where Hq is the Hubble constant at present. This can be rewritten as 



(a) = —[1 + zl)Dol / T^T-To- (3-16) 



-Dos Ho Jy^ D{z 

The right hand side of this equation depends on a only through Dqi^/ D[z). Following 
reasoning similar to that used in the proof of Eq. (p-5D , we obtain for zi^ < z < z^ 



^^(ai) < ^^("2) for ai < 02- (3-17) 
From Eqs. ( |34^ ) and ( FT7| ), we obtain 

— F; («i) < — f; («2) for ai < a2. (3-18) 

Therefore, the gravitational lensing event rate increases with a. 
(3) Di^s/DolDos 

Finally, we investigate the combination of distances appearing in the time delay. 



Dividing Eq. (3^4) by Eq. ( p-5D , we obtain 

— f; ("1) > — f; («2) for ai < 02- (3-19) 

Thus, the time delay decreases with a. 

As shown above, the three types of combinations of distances are monotonic func- 
tions of the dumpiness parameter. However, some of other combinations of distances 
are not monotonic functions of a, though these combinations may not be necessarily 
related with the observation. For instance, the combination -Dls / V cDqs / Hq is not 
a monotonic function of a. 
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3.4. Implications for cosmological tests 

We consider three types of the cosmological test which use combinations of 

distances appearing in gravitational lensing. Let us fix the density parameter in 
order to discuss constraints on the cosmological constant. 
(1) Dj^s/Dos 

The following relation holds 

^(Ai)<^(A2) for Ai<A2. (3-20) 
This is shown as follows. Let us define 

= p^^^. (3.22) 



and 



By the reciprocity, 11) we obtain 



Yxizi^) = i^X,(zL), (3-23) 

1 + 

which satisfies 

(f 3 

— ^n(^L) + + zi^fa^oY^izL) = 0. (3-24) 
dfL 2 

For Xi (i = 1, 2), the affine parameter Vi satisfies 



:i + zl)V^o^l(1 + ^l)'-A*^l(2 + zl) + (1 + zl)2. (3-25) 
We define the Wronskian as 

W(Y,,,Y,,) = Y,,'-X^-Y,/^. (3.26) 

Then, using Eq. ( 3-24| ), we obtain 

-^W{Yx„Yx,) <0 for Ai < A2. (3-27) 

Since Yx always vanishes at zi, = zs, we also obtain 

W{Yx,{zs),YxAzs)) = 0- (3-28) 
From Eqs. ( |3^ ) and (|]2|), we find 

W{Yx„Yx,)>0 for Ai<A2, (3-29) 

which can be rewritten as 

Ain^>0 for Ai<A2. (3-30) 

dzL 
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Since Y\ always becomes 1 + at the observer, we have 

n2(^L = o) 



In. 



0. 



(3-31) 



(3-32) 



ni(^L = o) 

Finally from Eqs. (|3-30D and (|3^ ), we obtain 

^ > 1 for Ai < As. 

Thus, Eq. ( F20|) is proved. 

Equations (|3-5| ) and ( ^-201 ) imply that, in a cosmological test using the bending 
angle, the cosmological constant estimated by use of the distance formula in the 
FLRW universe is always less than that etimated by use of the Dyer-Roeder distance 
(0 < a < 1). 
(2) Doi^Di^s/Dos 

Multiplying Eq. (|^) by 



we obtain 



^ol(Ai) < DolIAs) for Ai < A2, 



"(Ai) < — 7:; (A2) for Ai < A2. 



^os 



^os 



(3-33) 



(3-34) 



Equation (3-33) can be proved, for instance, in the foll owing manner: The Dyer- 
Roeder distance is written as the integral equation E^^'E^P 



D{z]a) = D{z;a = + 



i=l 



where Ki{y,z) is defined as 



Kiix,y) 



3_c_ 

2H'n 



dv 
dz 



{l-a)Q 



r dyKiiy,z)D{y;a = l), (3-35) 
Jo 



[l + x)^Dix,y;a = l) 



and 



Ki+i{x,y)= dzKi{x,z)Ki{z,y). 

J X 

From Eqs. (|3^ and (|3^ ), it is shown that for x < y 

Ki{x,y;Xi)<Ki{x,y;X2) for Ai < A2, 
where we have used the relation 

D{x, y;a = 1, Ai) < D{x, y;a = 1, A2) for Ai < A2, 



(3-36) 
(3-37) 

(3-38) 
(3-39) 



applicable in the FLRW universe. Using Eqs. (|3^ , (|3^ and (|3^ , and the 
positivity of Ki , we obtain Eq. ( |3-33 ) . 

From Eqs. ( |348|) and (|]3|), it is found that, in a cosmological test using the 
lensing events rate, the cosmological constant is always underestimated by use of the 
distance formula in the FLRW universe. 
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(3) Dj^s/DolDos 

When the time delay is measured and the lens object is observed, -Dql-Dos/^ls 



can be determined from Eq. (3^). On the other hand, when we denote the dimen- 
sionless distance between zi and Z2 as di2 = HqDi2/c, which does not depend on 
the Hubble constant, we obtain 



Then, Eq. ( ^-19 ) becomes 



(3-40) 



dohdoS f ^ ^ dohdoS ^ . r ^ /Q/|1\ 

— (ai) > — (02) lor ai < 02- (3-41) 

"LS "LS 



Thus, from Eqs .( 3^401 ) and ( 3-41 ), it is found that Hq estimated using the Dyer- 



Roeder distance decreases with a. Thus, the Hubble constant can be bounded from 
below when we have little knowledge on the dumpiness of the universe. The lower 
bound is given by use of the distance in the FLRW universe. On the other hand, 
since the combination -DlsZ-Dol-Cqs is not a monotonic function of the cosmological 
constant, the relation between the dumpiness of the universe and the cosmological 
constant is not simple. 

It should be noted that even the assumption of a spatially flat universe (J7 + 
A = 1) does not change the above implications for the three types of cosmological 
tests, since the cosmological constant affects the Dyer-Roeder distance formula only 



through the relation between z and v, Eq. ( 2-26 ). 

3.5. Evolution of dumpiness 

We have taken the dumpiness parameter a as a constant along the line of sight. 
However, as a reasonable extension of the DR distance, a can be considered as a 
function of the redshift in order to take account of the growth of inhomogeneities 
of the universe. In proving the monotonic properties, it has never been assumed 
that a is constant. Hence, all the monotonic properties and the implications for 
cosmological tests remain unchanged for the variable a{z). That is to say, when 
ai{z) < a2{z) is always satisfied iov < z < zs, all we must to do is to replace 



parameters ai and 02 with functions ai{z) and a2{z) in Eqs. (|3-5|), (3-18) and 



(3-19). In particular, when a{z) is always less than unity on the way from the source 
to the observer, both of the combinations of distances appearing in (1) and (2) are 
less, while the combination in (3) is larger than those in the FLRW universe. Then, 
the decrease in the bending angle and the lensing event rate, and the increase in the 
time delay hold even for a generalized DR distance with variable a{z). 
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§4. Average distances and the dispersions in inhomogeneous models 



In this section we describe the statistical behavior of angular diameter distances 
in inhomogeneous model universes at the stage of < z < zi(= 5). To derive the 
distances we use the light rays received by (or emitted backwards from) an observer at 
present (by solving the null-geodesic equation) in the universes which were produced 
numerically. 

4.1. Model universes and lens models 

We consider three background models with [Qq,Xq) = (1.0, 0), (0.2, 0.8) and 
(0.2, 0). They are denoted as S, L and O, respectively, which represent the standard 
model, a low-density flat model and an open model. The matter is assumed to 
contain particles consisting of galaxies and non- galactic clouds with equal mass m, 
but generally different sizes. The inhomogeneous models are given by the method 
of the A^-body simulation (using Suto's tree co de)El) in periodic boxes with particle 
number N = 32^. The initial particle distributions were derived using Bertschinger's 
software COSMICS^ under the condition that their perturbations are given as 
random fields with the spectrum of cold dark matter, their power n is 1, their 
normalization is specified as the dispersion erg = 0.94, and the Hubble constant is 
Hq = 100/iMpc~^ km s~^, where h = 0.5 for (1.0,0) and h = 0.7 for other models 
with Qq = 0.2. 

The box sizes for models S, L and O are 



Lo = a(io)^ = 32.5/i^\ 50/i-\ SO/i'^Mpc, 



and the particle masses are 



m{= pbqLq-^/N) = 2.90, 2.11, 2.11 x IO^/i^^Mq, 



(4-1) 



(4-2) 



respectively, where pBo is the background mass density, a{t) is the scale-factor, and 
/ is the comoving length. 

The particle size r<j (= a{t)Xs) is given in the form of softening radii, which have 
constant values when we calculate the gravitational potential for lensing. For we 
consider the following two (lens) models: 

Lens model i. All particles in the low-density models {Qq = 0.2) have = 
20/i~^kpc, 20% of the particles in the flat model (1.0, 0) have = 20/i^^kpc, and the 
remaining particles have r<j = 500/i~^kpc. Thus, practically, particles with = 0.2 
(which we call compact lens objects) play the role of lens objects. Their number 
density is much larger than the galactic density Qg ~ 0.02. 

Lens model 2. 10% of the particles in the low-density model (l7o = 0.2) and 2% 
of the particles in the flat model (1.0,0) have = 20/i^^kpc, while the remaining 
particles have = 500/i^^kpc. Therefore only galaxies corresponding to Qg = 0.02 
play significant roles as lens objects, and the remaining particles are regarded as 
diffuse clouds. 

The background line-element is 



ds^ 



-c^dt^ + a'^{t){dxf/ 



1 + k\{x) 
4 



(4-3) 
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and the Poisson equation and null-geodesic equation describing light rays are given 
in §2 of a separate paper (by Tomita, Premadi and Nakamura) of this volume. 

4.2. Angular diameter distances 

Here we treat the linear and area distances defined in §2 (Eqs. ( ^ ) and ( |2^ ). 
Let us consider a pair of rays received by the observer with the separation angle 9. 
By solving null-geodesic equations, the interval of the two rays at any epoch can 
be derived. If {Ax)± is the component of the deviation vector perpendicular to the 
central direction of the rays, the linear angular diameter distance -Dia is given as 



Aa = a{t)iAx)_ 



1 + ]k{x)^ 
4 



/0, 



(4-4) 



where the factor (1 — 2^) has been neglected, because <^ 1 locally. The above 
expression can be rewritten by use of (= qqx'^/Ro) as 



{l + z)F 



(4-5) 



1/3 



where F = 1 - i(i?o^o/c)2(l - - Ao)(y)^ ao = a{to) and Rq = Lq/N 

On the other hand, the area angular diameter distance DaA is given as follows 
using three rays (ray 1, ray 2 and ray 3) received by the observer, such that on the 
observer plane the two lines between ray 1 and ray 2, and between ray 1 and ray 
3 are orthogonal and have the same lengths (equal to the separation angle 9). If 
(Z\a;)^(i2), (/ia;)_|_(i3) and (Ax) ^(^23) are the components of the deviation vectors 
(between ray 1 and ray 2, between ray 1 and ray 3 and between ray 2 and ray 3 ) 
perpendicular to the central direction of the rays, we obtain 



D. 



aA 



a(t)[(Z\;r)^(i2) • (Z\a;)^(i3)]i/2 



l + ^Kixf 
4 



/9. 



(4-6) 



Using y* (= aQX^/Ro), this reduces to 



Rn 



{l + z)F 



-AyuAyisKAyu 



Ayi3 + Ay23){Ayi2 - Ayi3 - Ay: 



'23 j 



1/4 

/e, 

(4-7) 

where Ayp = |(Z\y)_|_(p)| with p = 12, 13, 23. 

In a previous paper (Tomita 0)) we investigated the behavior of Dia for the 
separation angle 9 = 0.005 — 20 arcsec in various model universes, and found the 
dependence of distances on 9 is small for 9 < 1.0 arcsec. Here we fix the separation 
angle to 9 = 1.0 arcsec and consider the difference between the linear and area 
distances and their dependence on the lens models 1 and 2. 

In the present lensing simulation we performed the ray-shooting of 500 ray bun- 
dles to derive Dia and DaA for each set of two lens models and three model universes. 
At the six epochs z = 0.5, 1, 2, 3, 4 and 5, we compared the calculated distances with 
the Dyer-Roeder distance and determined the corresponding value of the dumpiness 
parameter a as follows. In Ref. 17) we calculated a for < z < 5 in the above 
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three model universes and found that the angular diameter distance depends ap- 
proximately linearly on a (cf. Figs. 3 6 in Ref. 17). For |a — 1| ^ 1 linearity does 
not hold (cf. Eq. (3.35)), but most light rays arc in the neighborhood of a = 1, as 
is verified below. Hence we determined a from the calculated distance D\ (= Dia 
or -DaA) using the relation 



a 



{Da - Dbr)/{Dp - Ddr), 



(4.J 



where Ddr is the limiting Dyer-Roeder distance with a = 0, and Dp is the calculated 
Friedmann distance in the homogeneous case. This Dp is equal to the Dyer-Roeder 
distance with a = 1. Moreover, we consider the normalized distances defined by 



dA = Da/Df. 



(4-9) 




Fig. 1. The percentage {WON{a)/N) of the 
distribution of a in bins with the inter- 
val Aa = 0.4, for Dia in the lens model 1 
and model S with (i7o,Ao) = (1.0,0). Re- 
sults for z = 0.5, 1, 2, 3, 4 and 5 are denoted 
by dot-long dashed, dot-short dashed, long 
dashed, short dashed, dotted and solid 
lines, respectively. 



Fig. 2. The percentage {100N{a)/N) of the 
distribution of a in bins with the interval 
Aa = 0.4, for Dia in the lens model 2 and 
model S with (1.0,0). The lines have the 
same meaning as in Fig. 1. 



As a result of statistical analysis for this ray-shooting, we derived the average 
dumpiness parameter a, the average normalized distance dA, their dispersions {aa 
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Table I. The average dumpiness parameter a and its dispersion aa, and the average normalized 
distance dA and its dispersion ad for two lens models in model S with {Oo, Ao) = (1.0, 0). 









linear 








area 


{Da.A) 




lens 








di 




ol 


^ a 




a 




u.o 


1 HQ 

i.uy 




1 on 


U.Uio 


1 flQ 

i.uy 


U.oD 


l.UU 


U.UIO 




1 


1.03 


0.25 


1.00 


0.032 


1.04 


0.22 


1.00 


0.028 


1 


2 


1.02 


0.15 


0.99 


0.051 


1.03 


0.41 


0.99 


0.045 




3 


1.01 


0.11 


0.99 


0.063 


1.02 


0.10 


0.99 


0.056 




4 


1.01 


0.09 


0.99 


0.073 


1.02 


0.08 


0.99 


0.064 




5 


1.01 


0.08 


0.99 


0.080 


1.02 


0.07 


0.99 


0.071 




0.5 


1.09 


0.28 


1.00 


0.012 


1.07 


0.23 


1.00 


0.009 




1 


1.03 


0.14 


1.00 


0.019 


1.03 


0.13 


1.00 


0.017 


2 


2 


1.02 


0.09 


0.99 


0.029 


1.02 


0.08 


0.99 


0.027 




3 


1.01 


0.07 


0.99 


0.037 


1.02 


0.06 


0.99 


0.034 




4 


1.01 


0.06 


0.99 


0.043 


1.01 


0.05 


0.99 


0.040 




5 


1.01 


0.05 


0.99 


0.048 


1.01 


0.04 


0.99 


0.044 



Table II. The average dumpiness parameter a and its dispersion aa, and the average normalized 
distance dA and its dispersion ad for two lens models in model L with (i?o, Ao) — (0.2, 0.8). 









linear 


(Aa) 






area 


(/?aA) 




lens 


z 


a 


Oci 


dA 


ad 


a 




dA 


ad 




0.5 


1.09 


1.43 


1.00 


0.019 


1.02 


0.91 


1.00 


0.012 




1 


1.07 


1.15 


1.00 


0.057 


1.02 


0.69 


1.00 


0.034 


1 


2 


1.02 


0.68 


1.00 


0.114 


1.03 


0.49 


0.99 


0.083 




3 


1.01 


0.50 


1.00 


0.157 


1.03 


0.37 


0.99 


0.118 




4 


1.00 


0.39 


1.00 


0.188 


1.02 


0.29 


0.99 


0.138 




5 


1.00 


0.33 


1.00 


0.211 


1.02 


0.24 


0.99 


0.155 




0.5 


1.08 


0.43 


1.00 


0.006 


1.06 


0.29 


1.00 


0.004 




1 


1.04 


0.33 


1.00 


0.016 


1.03 


0.21 


1.00 


0.010 


2 


2 


1.01 


0.20 


1.00 


0.034 


1.01 


0.13 


1.00 


0.022 




3 


1.01 


0.15 


1.00 


0.047 


1.01 


0.10 


1.00 


0.032 




4 


1.01 


0.12 


1.00 


0.056 


1.01 


0.08 


1.00 


0.039 




5 


1.01 


0.10 


1.00 


0.064 


1.01 


0.07 


0.99 


0.045 



Table III. The average dumpiness parameter a and its dispersion a a, and the average normalized 
distance dA and its dispersion ad for two lens models in model O with (Jlo, Ao) = (0.2, 0). 









linear 


(Aa) 






area 


(^-aA) 




lens 


z 


a 


aa 


dA 




a 


fa 


dA 


fd 




0.5 


1.02 


1.99 


1.00 


0.020 


1.07 


1.24 


1.00 


0.012 




1 


1.07 


1.55 


1.00 


0.046 


1.12 


1.22 


1.00 


0.037 


1 


2 


1.06 


1.06 


1.00 


0.084 


1.09 


0.84 


0.99 


0.067 




3 


1.04 


0.89 


1.00 


0.113 


1.09 


0.74 


0.99 


0.094 




4 


1.03 


0.82 


0.99 


0.137 


1.08 


0.68 


0.99 


0.115 




5 


1.03 


0.75 


0.99 


0.156 


1.07 


0.62 


0.99 


0.129 




0.5 


1.08 


0.61 


1.00 


0.006 


1.04 


0.44 


1.00 


0.004 




1 


1.08 


0.49 


1.00 


0.014 


1.08 


0.41 


1.00 


0.012 


2 


2 


1.04 


0.31 


1.00 


0.025 


1.06 


0.26 


1.00 


0.021 




3 


1.03 


0.26 


1.00 


0.033 


1.06 


0.22 


0.99 


0.028 




4 


1.02 


0.24 


1.00 


0.040 


1.05 


0.20 


0.99 


0.034 




5 


1.02 


0.22 


1.00 


0.046 


1.05 


0.18 


0.99 


0.038 
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Fig. 3. The percentage (100A''(Q)/iV) of the 
distribution of a in bins with the interval 
Aa = 0.4, for DaA in the lens model 1 and 
model S with (1.0,0). The lines have the 
same meaning as in Fig. 1. 



Fig. 4. The percentage {WON{a)/N) of the 
distribution of a in bins with the interval 
Aa = 0.4, for DaA in the lens model 2 and 
model S with (1.0, 0). The lines have the 
same meaning as in Fig. 1. 



and cJd) and the distribution {N{a)) of a. Here, in order to study the frequency of ray 
pairs with a, we consider many bins with the interval Aa = 0.4 and the centers Oj = 
1.0 lb OAi [i = 0, 1, 2, ...). The number of ray pairs with — Aa/2 < a < + Aa/2 
is expressed as N{ai) and the total number of ray pairs is A^(= J2i ^i<^i))- In Tables 
I ~ III, {a,aa,dA, CTd) for and Z?aA in the two lens models are shown for models 
S, L and O, respectively. In Figs. 1 ~ 12, the percentages of the distribution of a, 
that is, 100A^(a)/A^ for D\a and DaA are shown for the above three models. The bar 
graphs in these figuresJiave the same meaning as the line graphs which were used 
in the previous paper. O In Figs. 1, 5 and 9, the distributions for six values of z 
are shown, and in the other figures those for only z = 1,3 and 5 are shown to avoid 
confusion. The following types of statistical behavior are found from these tables 
and figures : 

(1) In all cases, both the average values a and are nearly equal to 1, so that the 
average angular distance can be regarded as the Friedmann distance. 

(2) For each angular distance, the two kinds of dispersions have different behavior: 
cTq, increases in the order of O, L and S for the same value of z, and cJq, in a given 
universe model decreases with the increase of z. On the other hand, ad increases 
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Fig. 5. The percentage {WON{a)/N) of the 
distribution of a in bins with the interval 
Aa = 0.4, for Dia in the lens model 1 and 
model L with (0.2, 0.8). The Unes have the 
same meaning as in Fig. 1. 



Fig. 6. The percentage {WON{a)/N) of the 
distribution of a in bins with the interval 
Aa = 0.4, for Dia in the lens model 2 and 
model Ol with (0.2, 0.8). The lines have 
the same meaning as in Fig. 1. 



in the order of L, O and S, and a a in a given universe model increases with the 
increase of z. This behavior is connected with the situation that the change in the 
Dyer-Roeder distance corresponding to the change in a is small for a < 1. 

(3) Generally the dispersions for a-re larger than those for i^aA, and the ratios 
of two dispersions are ~ 1.2 for model S and ~ 1.6 for models L and O. These 
differences can be seen also by comparing Figs. 1 and 3, Figs. 2 and 4, and Figs. 
10 and 12. 

(4) The dispersions in the lens model 2 are smaller than those in the lens model 1. 
The ratios of the former to the latter are ~ 2/3, 1/3, 1/3 for universe models S, L 
and O, respectively. These differences can be seen similarly by comparing Figs. 1 
and 2, Figs. 3 and 4, and Figs. 11 and 12. In the two lens models of model S 
and in the lens model 2 of models L and O, the angular diameter distances can be 
regarded as the Friedmann distance, because of the small dispersions. In the lens 
model 1 of models L and O, however, we cannot always use the Friedmann distance, 
because of comparatively large dispersions. 

Using numerical ray-shooting in the A^-body-simulating clumpy cosmological 
models, we studied the statistical behavior of the angular diameter distances D\a 
and -DaA and determined the dumpiness parameter a by comparing it with the 
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Fig. 7. The percentage {WON{a)/N) of the 
distribution of a in bins with the interval 
Aa = 0.4, for DaA in the lens model 1 and 
model 02 with (0.2,0.8). The lines have 
the same meaning as in Fig. 1. 



Fig. 8. The percentage (100A'^(q)/A'^) of the 
distribution of a in bins with the interval 
Aa = 0.4, for DaA in the lens model 2 and 
model 02 with (0.2, 0.8). The hnes have 
the same meaning as in Fig. 1. 



Friedmann distance (a = 1) and the Dyer-Roeder distance (a > 0). Moreover, we 
studied the behavior of the normalized distance dx- The results show that all average 
values of a arc nearly equal to 1, the dispersions of the linear distance arc slightly 
larger than those of the area distance, and in the lens model 1 of models L and O, 
the dispersions are not so small, and we cannot use the Friedmann distance, while 
in the other cases we can use the Friedmann distance because of small dispersions. 

In the above averaging process, all light rays were taken into account. If we 
consider only weakly deflected light rays as contributing to weak Icnsing, the disper- 
sion (Jq will be slightly smaller than the values in the above tables. However, the 
contribution of strong lensing to aa is small because of its small frequency. 

Finally, we touch on the estimate of Icnsing correction to source magnitudes 
based on Tables I ~ III. Using aa in the area distance it is given by 

Amiz) = ^ log[l + 2ad\ = 2.18(7^ (4-10) 

for the source at z. Its values for the lens model (1,2) are 
Am{0.5) = (0.033, 0.020), (0.026, 0.009) and (0.026,0.009), 
Z\m(1.0) = (0.061, 0.037), (0.074, 0.022) and (0.081,0.026) 
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Fig. 9. The percentage (100A''(Q)/iV) of the 
distribution of a in bins with the interval 
Aa = 0.4, for Dia in the lens model 1 and 
model O with (0.2,0). The lines have the 
same meaning as in Fig. 1. 



Fig. 10. The percentage (100iV(a)/iV) of the 
distribution of a in bins with the interval 
Aa = 0.4, for Dia in the lens model 2 and 
model O with (0.2, 0). The lines have the 
same meaning as in Fig. 1. 



for models S, L and O, respectively. This Am (for the separation angle 6 = 1 arcsec) 
was derived independently of Am (for the separation angle 9 = 2 arcsec), shown 
in the §3 of Tomita, Premadi and Nakamura's paper in this Supplement, but their 
values at 2: = 1 are roughly consistent. The lensing correction to source magnitudes 
was also investigated by Holz in different lens models and inhomogeneous models. e2f 



§5. Shear effect on distances 

As we have shown in §2, the evolution of a cross sectional area of a light ray 
bundle is determined by Ricci and Weyl focusing along the trajectory of the ray 
bundle. Ricci focusing is a convergency effect due to matter in the ray bundle. On 
the other hand, Weyl focusing is a result of the tidal shear on the ray bundle induced 
by the inhomogeneous distribution of matter. 

One of the main difficulties in deriving a distance-redshift relation analytically 
for a realistic inhomogeneous universe lies in estimating the effect of Weyl focusing 
on ray bundles. 

Since the pioneering work of Gunn, cS* there has been a great deal of progress in 
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Fig. 11. The percentage (100iV(a)/iV) of the 
distribution of a in bins with the interval 
Aa = 0.4, for DaA in the lens model 1 and 
model O with (0.2,0). The lines have the 
same meaning as in Fig. 1. 



Fig. 12. The percentage (100iV(a)/iV) of the 
distribution of a in bins with the interval 
Aa = 0.4, for DaA in the lens model 2 and 
model L with (0.2, 0). The lines have the 
same meaning as in Fig. 1. 



constructing an analytical approach to investigate statistical quantities of realistic 
distances (e.g., dispersion and skewness of probability distribution of image magni- 
fications). Babul and Lee,E3 among others, examined lensing magnification effects 
on distances due to the large scale structures (> 0.5/i~^Mpc, where Hubble constant 
Hq = 100/ikm/sec/Mpc). They found that the dispersion in image magnifications 
is negligible even for sources at a redshift of 4. At the same time, they pointed out 
that the dispersion is very sensitive to the nature of the matter distribution on small 
scales. In their study, the effects of Weyl focusing were neglected based, on a numer- 
ical study by Jaroszyhski et al.cS' in which the lensing magnification effects due to 
large scale structure (> l/i~^Mpc) in cold dark matter models were examined, and it 
was concluded that Weyl focusing has no significant effect on image magnifications. 
FriemanliZP improved Babul and Lee's study to reflect the recent developments in 
numerical and observational studies of large scale structure. Although he took small 
scale nonlinear structure into account, Weyl focusing effects were neglected without 
any reasonablcLljasis. 

NakamuraS examined the effects of shear on image magnification in the cold 
dark matter model universe with a linear density perturbation. He found that the 
effect is sufficiently small and concluded that Weyl focusing can be safely neglected 
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for a light ray passing through linear density inhomogeneities. 

The above cited studies mainly focus on large scale inhomogeneities, whereas the 
effects of small scale objects, such as galaxies and clusters of galaxies, have not been 
satisfactorily taken into account. It is, however, not clear whether the Weyl focusing 
effect due to small scale inhomogeneities has a significant effect on distances. In 
this section, we discuss theWeyl focusing effect due to small scale inhomogeneities, 
mainly following Hamana. C§' 

5.1. Basic equations 

We first derive^an evolution equation of lensing magnification frorn thfi null 
geodesic equation, e3) which is equivalent to the optical scalar equations, Ea) 'ESP and 
is convenient to examine gravitational lensing effects. 

We rewrite the cosmological Newtonian metric ( |1-2D as 

ds^ = a^{rj) + 2<P)drf' + (1 - 2'!>)-iijdx' dx^] , (5-1) 

where ry is a conformal time: drj = cdt. We write the above metric as g^^y = a'^g^i^. 
Since the light cone structure is invariant under the conformal transformation of the 
metric, in the following we work in conformally related g world. 

Let us consider an infinitesimal bundle of light rays intersecting at the observer. 
We denote a connecting vector which connects the fiducial light ray 7 to one of its 
neighbors as All gravitational focusing and shearing effects on the infinitesimal 
light ray bundle are described by the geodesic deviation equation, 

^ = -R^au^^k^", (5-2) 

where k" = dx^/dX, and A is the affine parameter along the fiducial light ray 7. 
We introduce a dyad basis e^"^ {A, B, C, ... = 1,2) in the two-dimensional screen 
orthogonal to k^ and parallel-propagated along 7. The screen components of the 
connection vector are given by 

= e/C^ (5-3) 

From the geodesic deviation equation ( ^•2[ ), one can immediately find that Y"^ sat- 
isfies the Jacobi differential equation 

^g^^r.«r«, (5.4) 

where Tab is the so-called optically tidal matrix. From the metric ( |5-1| ), up to 
first order in 'P, this matrix is given by 

where I is the 2x2 identity matrix, and TZ and represent the Ricci and Weyl 
focusing induced by the density inhomogeneities, respectively: 

7^ = Z\(% = ^^^, (5-6) 
2 a 

^=^11-^22 + 2^^,12. (5-7) 
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Here A^^^ is the Laplacian operator in the spatial section, and 6 is the density 
contrast defined hy 6 = {p — pb)/pb, where pb is the mean matter density. Owing to 
the linearity of (|5-4D , the solution of can be written in terms of its initial value 
dY^/dX\x=Q = "d^ and the A-dependent linear transformation matrix T^ab can be 
written as 

y^(A) =P^b(A)^?^. (5-8) 
Substituting the last equation into the Jacobi differential equation , we obtain 



Now, we derive an evolution equation of the lensing magnification matrix relative 
to the smooth Friedmann distance from ( ^-91) . First, we write ( |5-5D as T=T^^^+ST, 
with T^^^ = —KI, and ST is the second term in { f>-5\ ). In the homogeneous case, ST 
is vanishing, and the solution of X> is VabW = Df{X)6AB, where Df is, of course, 
the standard angular diameter distance in the background Friedmann universe. 

It is natural to define the lensing magnification matrix relative to the corre- 
sponding Friedmann universe as 

Differentiating Mab twice with respect to A and using (p-9\), one finds 



d'^MAB 2 dPf dMAB , 

= -D}liX^Or + ^^^^-^^^^ (5-11) 

With the initial conditions M(A)|a=o =^ and dM {\)/d\\x=o =C>,0 the last 
equation can be written in the integral form 

Mab{\) = 6ab + r dX' ^^^^'^'l^^^^'^ 6Tac{X')Mcb{X')- (5-12) 
Jo 

This is the general form of the evolution equation of the lensing magnificatioiLmatrix 
relative to the Friedmann distance in multiple gravitational lensing theory. EP Note, 
in general, this equation is not an explicit equation for Mab^ since it involves an 
integration over the optical tidal matrix evaluated on the light ray path, such that one 
first has to solve a null geodesic equation. Sjjice, for almost all cases of cosmological 
interest, the deflection angle is very small, ^ we will neglect the deflection of light 
rays. 

5.2. Order- of -magnitude estimate 

We now examine the magnitude of lensing effects on light ray bundles due to ran- 
domly distributed virialized objects (e.g. galaxies and clusters of galaxies) adopting 
an order-of- magnitude estimate. eP'lIP As can be seen in Eq. ( |5-12D , if the magnitude 
of the components of the matrix, J dX'Df{X — X')Df{X')/Df{X)6'T{X') is small, the 
magnitude of the lensing effects is dominated by these terms. We, therefore, examine 
the magnitude of these terms. For simplicity, we denote these terms as 5Ai. 
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Supposing that lensing objects are randomly distributed and that each has a 
mass M = 2av'^l/G, where is the one-dimensional velocity dispersion of the lens 
objects, and I is a characteristic comoving size of a lens object. Hence the mean 
comoving number density of the lens objects is 



ZQlHo^ 1 1 



(5-13) 



where J?^ is the density parameter of lens objects defined by J7l = pl/ {"^Hq /SttG), 
with PL is the mean density of the lens objects. Thus, the mean comoving separation 
distance is 

^ ^^"^ .aH\\ (5.14) 



Then for a geodesic affine comoving distance of A, the light ray gravitationally en- 
counters such objects Ng = X/vq times on average. At each encounter, the contribu- 
tion to the lensing magnification matrix is 



AM 



47r 



^rp Df{Xd)Df{Xs-Xd) 

52 Dj{Xs) 
2 1 Df{Xd)Df{Xs-Xd) 

ro DfiXs) 



(5-15) 



where Df{Xi — Xj) is the comoving angular diameter distance, the subscripts d and s 
indicate the lens and source, respectively, and b is the comoving impact parameter. In 
the above expression, we have assumed that the mean comoving impact parameter is 
of order tq. Since the sign of each contribution will be random, the total contribution 
to the lensing magnification matrix is 

SM ~ AMx/K 



7V^ 
4 



'HoDf{Xd)Df{Xs-Xd) 



\ c Df{Xs) I 



c 



(5-16) 



where i/ is a compactness parameter of a lens object defined by = 47r(c7^/c)2/~^ 
X {c/Hq). The contribution from direct encounters can be similarly estimated by 
noting that the average number of encounters is 



fx 3f2LHo^ 1 

^^d — — = — 77^ ^t^) 



'0 



IGtt ct^^ 



with each encounter contributing 



AMd = 4Tr{^] - 



ay\^ lDf {Xd)Df{Xs-Xs 
DfXs 



(5-17) 



(5-18) 
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with random sign. The result turns out to be the same as that of gravitational distant 
encounters, given by Eq. (5-16). The comoving affine distance A becomes c/Hq at 
the source redshift Zg ~ 3, and the averaged value of the distance combination over 
the lens redshifts is of order 



Accordingly, we find that the magnitude of the total contribution of the lensing 
effects to the lensing magnification matrix scales as ~ O.lV ^L\/v for the source 
redshift Zg ~ 3. We have the relation 17/, < i7o by definition, and l7o appears to be 
less than unity. Thus i?/, < 1. On the other hand, v <\ for galaxies and clusters 
of galaxies. Therefore a typical value of gravitational lensing effects on the lensing 
matrix can be expected to be O(O.l) or smaller for a majority of random lines of 
sight. 

The lensing magnification factor of a point like image is defined by the deter- 
minant of the lensing magnification matrix. Taking the determinant of the lensing 
magnification matrix ( [5T2D , and expanding it in powers of 5A^, one can easily find 
that the leading term of Weyl focusing effects is of order bM?' . On the other hand, 
that of the Ricci focusing term is of order bM.. Since we have seen that a typical 
value of bM. is expected to be 0(0.1) or smaller, we can conclude that, at least 
from a statistical point of view, Weyl focusing has no significant effects on the image 
magnifications or equivalently on the distances. 

5.3. Numerical investigation 

The above argument may sound too naive. One of the authors (T.H.), numeri- 
cally investigated Weyl focusing effects on image magnifications by using the multiple 
gravitational lens theory. e5P He focused on gravitational lensing effects due to small 
scale virialized objects, such as galaxies and clusters of galaxies. He considered a 
simple model of an inhomogeneous universe. The matter distribution in the universe 
was modeled by randomly distributed isothermal objects. He found that, for the 
majority of the random lines of sight, Weyl focusing has no significant effect, and 
the image magnification of a point like source within a redshift of 5 is dominated by 
Ricci focusing. He also found that his result agrees well with the order-of-magnitude 
estimate given above. 

To summarize, we conclude that, except for a statistically very rare kind of light 
ray, Weyl focusing has no significant effect on image magnifications or equivalently 
on the distances. 



Lensing observation in inhomogeneous universes was discussed in §3, based on 
the so-called Dyer-Roeder distance, in which one of the main assumptions is neglect- 
ing Weyl focusing. Such a neglection seems correct in our universe, as was shown in 
the §5. 

The average angular distances in inhomogeneous model universes are the Pried- 




(5-19) 



§6. Concluding remarks 
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mann distances, as was shown in §4, but individual ray bundles have various values 
of dumpiness parameters a because of their dispersions. The observational quanti- 
ties are sensitively dependent on a, as was shown in the §3, and so they may have 
dispersions similar to a. 

The difference between linear and area angular diameter distances, which is 
caused by the shear, is generally small in accord with the result in §5, even though we 
considered small-scale inhomogcncitics, but the difference in the low-density models 
is found to be larger than that in the Einstein-de Sitter model. This implies that 
the shear effect is comparatively larger in the low-density models. 
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